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1 Abstract. Following our previous work on metal abundances of a double damped Lya system with a line-of-sight 

separation ~ 2 000 km s -1 (Ellison & Lopez 2001), we present VLF UVES abundances of 3 new systems spanning 
a total of ~ 6 000 km s _1 at z ~ 2.5 toward the southern QSO CTQ247. Fhese abundances are supplemented with 
echelle observations of another 'double' damped Lya system in the literature. We propose a definition in terms of 
velocity shift of the sub-class 'multiple damped Lya system', which is motivated by its possible connection with 
large-scale structure. We find that the abundance ratio [S/Fe] is systematically low in multiple systems compared 
with single systems, and with a small scatter. The same behavior is found in 2 more single DLA systems taken 
from the literature that show evidence of belonging to a galaxy group. Although [Si/Fe] ratios are also generally 
lower in multiple DLAs than in single DLAs, the effect is less striking since the scatter is larger and there are 
' a number of low [Si/Fe] DLAs in the literature. We suggest that this can be explained with a combination of 

detection bias and, to a lesser extent, the scatter in ionization corrections for different absorbers. We investigate 
whether consistently low a/Fe ratios could be due to dust depletion or ionization corrections and find that the 
former effect would emphasize the observed trend of low a/Fe in multiple systems even further. Ionization may 
have a minor effect in some cases, but at a level that would not change our conclusions. We thus conclude that 
the low a/Fe ratios in multiple DLAs have a nucleosynthetic origin and suggest that they could be explained by 
reduced star formation in multiple damped Lya systems, possibly due to environmental effects. There seems to 
be independent evidence for this scenario from the mild odd-even effect and from the relatively high N/a ratios 
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we observe in these multiple systems. 
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1. Introduction dust content, photoionization and metallicity evolution. 

With the present sample of some 60 abundance mea- 
Dampcd Lyman Alpha systems (DLAs) m the spec- , uj+ju u -a a t 

. . . .„ v ' . surements some broad trends have been identified (e.g., 

tra of high-redshift quasars have promised to be a | Prodlaska fc Wolfc 2002) . but the global picture of star 

powerful probe of early metal enrichment in young formation (SF ) history deduced from these patterns of el- 
galaxies. However, despite the considerable effort in- e mental abundances -likely affected by dust depletion pat- 
vested to develop this technique, its full potential has terng j H ou, Boissier k Prantzos 2001 ^ - is still unclear. In 

yet to be realized. Although extensive observing cam- . , — , ~~r- — — r~ ■ ,£,„,. , . 

J iii nil particular, extracting coherent scenarios of SI histories in 

paigns at intermediate resolution have swollen the list . ,. . , , ^ A , , c , , . , , 

r ° . . ... ., , individual DLAs, or subsets of absorbers remains a cnal- 

of known DLAs (Elli son et al. 2001bl ILopez et al. 20011 , , , , , 

' ' ' — — -I lengmg prospect with current data. 



and IPeroux et al. 20011 are some of the most recent DLA . . ,. , , . ,. , . , . , , , . .. 

1 — ■ ■ ■ 1 One of the great advantages of studying high rcdshitt 

surveys), only a fraction have been followed- up at high ,. . , ,. . , , £ , ,. ,. 

, / ,.i ., • • . galaxies m absorption is the lack of selection bias asso- 

resolution to study their chemical properties, investigate . , i .,i i • i j- i t ± j 

J f t- i b ciated with cosmological dimming: galaxies are selected 

Send offprint requests to: S. Lopez based onl y on S as cross section, irrespective of their intrin- 

* The work presented here is based on data obtained at the sic luminosity. However, this feature may turn out to be a 

ESO Very Large Telescopes, programs 68.A-0361(A) and 66. A- double-edged sword in that extracting coherent chemical 

0660(A). enrichment patterns is almost certainly hindered by the 
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Table 1. VLT Observations of CTQ 247 



Mode 


Wavelength 


Exp. Time 


Observing 


Date 




M 


[sec] 






FORS2 (GRISM 600B) 


345-590 


600 


Dec. 18 


2000 


FORS2 (GRISM 600R) 


525-745 


400 


Dec. 18 


2000 


UVES Dichroic (390+580) 


328-445,476-684 


13 500 


Oct. 10, 11 and Nov. 17, 


2001 


UVES Dichroic (437+860) 


376-500,667-1040 


10 800 


Nov. 17 


2001 



mixed morphological representation. Undoubtedly, some 
of the scatter in relative abundances comes from vari- 
able dust depletion, although system-to-system differences 
persist even after correcting for grain fraction (Vladilo 
2002). In order to compare elemental trends in DLAs 
in a similar manner to, for example, local stellar abun- 
dance studies, a significant refinement in our approach is 
required. Such a refinement might include isolating par- 
ticular subsets of DLAs based either on their morphology 
or some other observed property (e.g. spin temperature; 
IChengalur fc Kanekar 2000| ). 

In this paper, we investigate whether a distinct enrich- 
ment trend exists for absorption line systems with close 
companions in velocity space, which we will refer to as 
'Multiple Damped Lyman Alpha systems', or MDLAs. 

2. The Sub-class 'MDLA' 

Lopez et al. (2001) identified the existence of three 
damped Lya absorbers spanning a total of v = 



5 900 km s" 



at z 



2.5 toward the southern quasar 



CTQ 247, the first ever example of a triple DLA; Ellison 
et al. (2001b) discovered a double absorption systems for 
which v = 1800 km s" 1 at z ~ 2 toward Q2314-409; 
a similar velocity span as in CTQ 247 is covered by the 
two DLAs toward Q2359-02 at z ~ 2 l|Wolfe et al. 1986|) . 
Statistically, these are extremely improbable events if the 
proximity of the systems is purely by chance, and only 
these few cases are reported in the literature. It is there- 
fore feasible to postulate that in these cases the multiple 
absorption systems may be in some way related, rather 
than chance alignments. However, such large velocities 
are difficult to reconcile with today's galaxy groups and 
clusters -i.e., virialized entities- or even with superclus- 
ters. However, at z > current work on small angular 
fields shows evidence for very large structures on comoving 
scales as large as d ~ 100 h^Q Mpc (Av ~ 10 000 km s _1 if 
purely due to Hubble flow) at various redshifts. Evidence 
for large structure on many Mpc scales has been found 
in absorbing gas both across (e.g., |Williger et al. 2 002 
d ~ 100 h^ 1 Mpc at z ~ 1.3 ) and along the line of sight 
HQuashnock, vanden Berk fc York 1996) d ~ 100 hr 1 
Mpc at z ~ 3), in super-clusters (d ~ 20 hr 1 Mpc 
at z ~ 0.8; Hai nes, Campusano fc Clowes 2003| l and in 
Lyman-break galaxies (Steidel et al. 1998; d ~ 10 
Mpc at z ~ 3 ). Large structure has also been reproduced 
in CDM N-body simulations (d ~ 20 /i" 1 Mpc at z ~ 1; 



lEvrard et al. 2002)) . It therefore seems a feasible (but not 
exclusive) possibility that absorption systems separated 
by several thousand km s _1 may be associated with very 
large scale structure. 

Motivated both by the uniqueness of MDLAs and the 
possibility of large-structure, we have embarked on a pro- 
gram to measure abundances in these systems. We will 
define an MDLA as 2 or more DLAs spanning a velocity 
range 



500 < Av < 10 000 km s~ 



(1) 



While the upper limit is set to match the limits of the 
largest structures known at high redshift, the lower veloc- 
ity roughly corresponds to the largest asymptotic veloci- 
ties in massive spirals (e.g. |Rigopoulou et al. 2002) ) 1 . Note 
that this a working definition only; although it might be 
possible that MDLAs probe some kind of large structure 
(but see §[5J, our definition does not necessarily imply that 
they host the strong C iv absorption used by Quashnock, 
vanden Berk & York (1996) or they are associated with 
the LBGs identified by |Steidel et al. (1996) 

From an instrumental point of view, identifying 
MDLAs in Lya -in contrast to using metal lines- will 
certainly limit the identification of systems separated by 
Av 1 000 km s _1 . If the discovery spectrum has poor 
S/N, obtaining line parameters for close MDLAs will re- 
quire both the detection of metal lines, and a good con- 
tinuum estimation 2 . 

For AT (Hi) criterion, we will relax the canonical limit 
to include systems with A 7 " (Hi) > 10 20 cm~ 2 . Not only 
do systems down to this limit continue to exhibit clear 
damping wings, but in the following we show that in such 
range ionization does not yet play a significant role when 
deriving abundances from low-ions. 

The first in-depth study of an MDLA was presented by 
Ellison & Lopez (2001; hereafter Paper I). There, the two 
DLAs toward Q2314-409 were noted to have low [S/Fe] 
ratios, a result which we suggested was due to the pos- 
sible impact of environment on the chemical abundances 

1 This lower limit will probably exclude the majority 
of bound galaxy satellites since most local spirals' satel- 
lites have velocity differences of less than 300 km s _1 , 
|Zaritsky et al. (1997)| 

2 Two AT (Hi) = 2 x 10 20 cm -2 DLAs separated by An = 
1000 km s^ 1 will mimic a single system in FWHM ~ 5 A data 



if S/N £ 20, but for Av 
S/N > 50. 



500 km s they will do so unless 
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in DLA protogalaxies. In the present work, we extend the 
study of abundances in MDLAs by investigating the triple 
DLA toward CTQ247 with new echelle data, and includ- 
ing literature abundances of the double system toward 
Q2359— 02. These new abundances reinforce our previous 
suggestion of peculiar relative abundances in MDLAs. 

After presenting the new data in § [3J the discussion 
of the abundances is given in § where we attempt a 
statistical comparison between the 2 populations, and as- 
sess possible systematic effects that may bias our result. 
A discussion on the nature of MDLAs is outlined in § 

3. Data Analysis 

3.1. Observations and Data Reduction 

CTQ 247 was observed in service mode with the UVES 
instrument at the ESO Kueyen telescope in October and 
November 2001 under good seeing conditions (0.6" — 1.0"). 
With the 390+580 and 437+860 dichroic modes we cov- 
ered from 328 to 1020 nm with two gaps at 576-583 nm 
and 852-866 nm. The exposure times were 13 500 seconds 
for dichroic 1 and 10 800 seconds for dichroic 2 (Table QJ. 

After the usual fashion of bias-subtracting and flat- 
fielding of the individual CCD frames, the echelle orders 
were extracted and reduced interactively with the UVES 
pipeline routines ( Ballcstcr et al. 2000J1 . The reference Th- 
Ar spectra used for wavelength calibration were taken af- 
ter each science exposure. The wavelength values were 
converted to vacuum heliocentric values and each order 
of a given instrumental configuration was binned onto a 
common linear wavelength scale of 0.04 A pixel -1 . The 
reduced orders were then added with a weight according 
to the inverse of the flux variances. Finally, the flux values 
were normalized by a continuum that was defined using 
cubic splines over featureless spectral regions. The spec- 
tral resolution is FWHM ~ 6.7 km s _1 , while the typical 
signal-to-noise ratio per pixel is S/N~ 35 — 40. 

In addition to the UVES data, a lower resolution 
(FWHM w 3 A) spectrum of CTQ 247 was obtained 
on December 18 2000 using FORS2 at the ESO Kueyen 
Telescope. This spectrum was used to better define the 
quasar continuum in the spectral region around the 
damped Lya lines. 

3.2. Column densities 

We used FITLYMAN and VPFIT 3 to fit the line profiles 
with theoretical Voigt profiles. All fits were unconstrained 
in redshift, Doppler width and column density, unless oth- 
erwise stated. In general, we prefer this approach over 
the apparent optical depth method (AODM; Savage & 
Scmbach 1991) when velocity components are not resolved 
as is sometimes the case. In addition, when dealing with 
transitions that lie in the Lya forest, such as the Sn 
triplet, fitting provides important information to the ex- 
tent of possible blending. Despite this, however, we did use 

3 Available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 



the AODM for transitions where the fit failed due to line 
blending of many components or poor S/N. We adopted 



the up-to-date /-values listed in Prochaska et al. (2001) 



0. 

□ 



and the solar abundances by [Grevesse fc SauvaT (1998) 
with updates for O and N by |Holweger (2001)| Table 
lists the column densities and derived abundances of all 
elements covered by our observations. 

The Hi column density was determined by fitting 3 
components to the damped Lya and Ly/3 profiles in the 
FORS2 spectrum. The choice of the low resolution spec- 
trum reduces the uncertainties inherent to determining 
iV(Hl) from echelle data, such as badly defined contin- 
uum where the damping wings extend over more than one 
order. We obtained iV(Hi) = 10 21 13 , 10 21 09 , and 10 20 47 
cm -2 , for the 3 systems centered at z a b s = 2.5505, 2.5950 
and 2.6215. The theoretical profiles are superposed on the 
data in Fig.^ The internal fit errors were 0.02 dex for all 
3 H I measurements but we believe a more realistic error 
is given by ±0.1, which is shown by the dotted curve in 
the figure. Although the noisy data around Ly/3 constrains 
Hi only marginally, this error seems quite safe, given the 
higher S/N at Lya. 

We next give a short description of the metal line fits 
in each of the 3 DLAs which we henceforth refer to as 
CTQ 247A,B and C. 



CTQ247A (z abs = 2.5505): For this DLA we have identi- 
fied nonsaturated transitions of Nl, Sin, Sn, Cm, Fell, 
and Znii, all of which show 4 velocity components. More 
(weaker) components are revealed by stronger transitions, 
e.g. Fen A2344, but these do not contribute significantly 
to the total column densities. Therefore, in order to avoid 
introducing too many parameters at intrinsically weak 
transitions in relatively noisy parts of the spectrum, we 
proceeded by fitting 4 components. Fig. |21 shows the non- 
saturated transitions used and their fits. From an analysis 
of the weaker components we estimate that the obtained 
column densities encompass ~ 90 % of the total column 
density; exclusion of weaker components will lead to only 
a ~ 0.05 dex difference at log N = 15 cm -2 . The Sn 1260 
triplet lines fall on the red wings of the damped Lya lines; 
consequently, they were fitted simultaneously with the H I 
fixed at the values found for the FORS data, hence the 
sub-unity values of the continuum in Fig. [5] In addition, 
we determine an AODM measurement for Mnii and an 
upper limit for P II. 

CTQ247B (z abs = 2.5950): This DLA has a relatively sim- 
ple velocity structure. Although the strongest Sin and 
Fe II transitions have 5 velocity components, these span 
only ~ 80 km s _1 . As expected for such a high Hi, Ol is 
saturated and we consider it no further here. Conversely, 
N I which is often barely detected exhibits a large equiva- 
lent width in this system. In fact, two of the A1134 triplet 
transitions are saturated. However, the triplet is relatively 
unblended, as are two of the S n triplet lines. Fig. shows 
the nonsaturated transitions used and their fits. There is 
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Fig. 1. Normalized flux and la VLT FORS2 spectrum of CTQ 247 at FWHM = 3 A resolution showing the Lya and 
Ly/3 lines of the triple DLA. Overlaid is a synthetic Voigt profile of three systems with 7V(H 1)= 10 2113 , 10 21 09 , and 
10 20.47 cm -2 ; cen tered at z abs = 2.5505, 2.5950, and 2.6215. The dotted lines are ±0.1 dex deviations for all AT(Hi) 
values. 



an absorption feature at the expected position of P II A963 
but since we suspect contamination by a Lya forest line we 
only provide an upper limit on N(P n). As for CTQ247A, 
we also provide an AODM measurement for Mnii. 

CTQ247C {z abs = 2.6215): This is the lowest Hi DLA 
of the triplet. For this DLA the bulk of the low ions are 
divided into 2 'sub-clumps' of 3 and 2 velocity components 
separated by ~ 120 km s _1 . Fig.0]shows the nonsaturated 
transitions used and their fits. Since this is the DLA with 
the lowest H I of the 3, fewer species are detected and limits 
are not very stringent. For example, for sulphur our 3cr 
detection limit is log N(S n) < 14.3. We have attempted 
a fit for the Ni triplet, but a combination of many weak 
components and Lya blending makes this a very uncertain 
result and we regard it as an upper limit due to possible 
contamination. On the other hand, the low Hi allows us 
to profile fit 2 O I transitions, for which each one of the 
5 components shows at least one nonsaturated line. Also, 
we report one of the few existing measurements of Al II in 
DLAs (line positions were held fixed at the Sill values), 
and one AODM measurement for Al in. 



4. Abundances in MDLAs 

In Paper I we focused on the relative a to Fe peak abun- 
dances in order to establish whether MDLAs may ex- 
hibit distinct SF histories. Similarly, in Fig. the solid 
circles depict [a/Fe] vs. [Fe/H] for the 3 MDLAs pre- 
sented in this paper (CTQ247A, B, and C), the dou- 
ble DLA_Q2^14-409A 2 _13(Paper I), and Q2359-02 A, 
B UProchaska fe Wolfe 1999)l . The stars are DLAs taken 
from the literature. Due to the potential problem of blend- 
ing between S lines and the Lya forest, we have compiled 



only values obtained from echelle data. In cases where 
numbers are present in several references, we have pre- 



ferred the values by Prochaska et al. (2001) for consis- 
tency. 

In the following we show statistically that the sub- 
class MDLA is drawn from a different distribution to sin- 
gle DLAs. We then discuss to which extent dust and ion- 
ization effects might produce such distinction, and what 
possible systematic effects may be present in the 2 sam- 
ples. 

4.1. Low a/Fe in MDLAs 

As expected in DLAs, all points in Fig. fall over the 
[a/Fe] = mark, as both S and Si are a-chain elements 
which are created mainly in massive stars on much shorter 
scales than the long-lived iron producing stars. The ra- 
tio of [a/Fe] is usually interpreted as an indicator of 
star formation history and the rate at which the differ- 
ent elements are released into the ISM. In addition to 
the MDLAs in Fig. we have also plotted DLAs with 
other neighboring galaxies (usually detected via Lyman 
break imaging) in the field. These systems are Q0201+112 
(Elli son et al. 2001a|l which has 4 identified galaxies in 
the vicinity, and Q00 00-262 l|Molaro et al. 200 1)1 for which 
Steidel et al. (1996) found 2 Lyman break galaxies within 
Az = 0.05 of the absorber. 

The remarkable property to note in Fig. 03 is that all 
MDLAs systematically show low [S /Fe] values with a small 
scatter (< 0.2 dex). To quantify this impression, we have 
performed a modified KS test that provides the likelihood 
of MDLAs with [S/Fe] measurements to be drawn from 
a different distribution to single DLAs. We created 1000 
realizations of a simulation which recreates the [S/Fe] dis- 
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CTQ247A CTQ247B CTQ247C 



X 


NfX] 


[X/H1 


NfX) 


IX/H1 


NfX) 


IX/H1 


ti I 


o i i o _i_ nin 
21. lo ± U.1U 




oi nn _i_ nin 
21.U9 ± U.1U 




on /i T _i_ nin 

20.47 ± 0.10 




N I 


14.55 ± 0.03 


n r* i i a 1 ri 

— 2.51 ± 0.10 


15.07 ± 0.02 


i n f i a 1 a 

— 1.95 ± 0.10 


< 14.36 


< —2.04 


U I 










i r i n l n nn 

15.19 ± U.U2 


o no _i_ n 1 n 
— 2.1)2 zb U.1U 


Alii 










12.88 ± 0.07 


— 1.97 ± 0.12 


Aim 


12.81 ± 0.05° 




12.63 ± 0.04 a 




12.33 ± 0.03° 




Sin 


15.32 ± 0.04 


— 1.37 ± O.ll 


15.59 ± 0.03 


— 1.06 ± 0.10 


i * k r\r\ i /"\ / \ 

13.99 ± 0.06 


— 2.04 ± 0.12 


P ii 


< 13.15 


< — 1.54 


< 12.98 


< —1.67 






S ii 


14.82 ± 0.06 


— l.5l ± 0.12 


15.19 ± 0.05 


— 1.10 ± 0.11 


< 14.34 


< —1.33 


Cm 


13.20 ± 0.03 


-1.62 ± 0.10 


13.37 ± 0.04 


-1.41 ± 0.11 






Mnll 


12.86 ±0.02 a 


-1.80 ±0.10 


12.87 ±0.01 a 


-1.75 ±0.10 






Fen 


14.95 ±0.06 


-1.68 ± 0.12 


15.15 ±0.02 


-1.44 ±0.10 


13.60 ±0.02 


-2.37 ±0.10 


Fein 






< 13.66 








Nin 


13.94 ±0.02 


-1.44 ± 0.10 


13.86 ±0.31 


-1.48 ±0.33 






Znn 


12.44 ± 0.05 


-1.36 ±0.11 


12.68 ±0.02 


-1.08 ±0.10 







Table 2. Abundance measurements (and 3cr upper limits) for CTQ247. a uses AODM; b N(Al) = AT(Aln) + AT(Alm) 



tribution of literature DLAs and MDLAs including the 1 
a quoted errors. For each DLA, an error was drawn at 
random from a Gaussian distribution and added to the 
observed value 4 . A KS test was then run on each of the 
1000 datasets and the likelihood that the KS probability 
Pks was greater than a certain confidence level was calcu- 
lated from the number of realizations. In this fashion, we 
overcome 2 problems associated with a normal KS statis- 
tics: the small number of MDLAs we have at our disposal 
and the exclusion of measurement errors. 

According to our simulations, the likelihood p that 
the two samples are drawn from different distributions is 
p = 1.00 for Pks > 95 % (i.e. every one of the 1000 realiza- 
tions gives a P K s > 95 %), p = 1.00 for P KS > 98 %, and 
even as high as 0.97 for Pks = 99 %. This indicates that 
the consistently low [S/Fe] values in MDLAs and other 
multiple systems are statistically significant. Since Si and 
S are both produced in oxygen burning reactions, [S/Si] 
= in all observed Galactic disk stars IjChen et al. 2002)l 
and we may expect to see similar patterns between S and 
Si in DLAs. However, the MDLA [Si/Fe] values have a 
slightly larger dispersion and there are a number of low 
literature [Si/Fe] values, which renders the two popula- 
tions less distinguishable. This may be explained because 
low iV(Sn) lines are easily contaminated in the Lya for- 
est, whereas low A^Sin) can be much more easily mea- 
sured. If we consider only the literature DLAs with Si and 
S measurements, the KS probability that the two [Si/Fe] 
samples are consistent is 11%, as opposed to 32% if all 
are included. There is an obvious outlier in the [Si/Fe] 
distribution associated with system Q2359— 02A; in the 
following section we argue that this system has anoma- 
lously large dust content compared with the rest of the 
MDLAs and literature DLAs. Discounting this exception- 



4 We treat the upper limit toward DLA B in Q2314-409 
as a detection; a lower value than this does not affect the KS 
probability 



ally dust depleted system reduces the KS probability even 
further of all [Si/Fe] measurements to 9%. Therefore, al- 
though there is a hint that [Si/Fe] may be low in MDLAs, 
the evidence is not as convincing as that for S. We explore 
this further in the following section. 

4.2. How dust and/or ionization might affect the 
observed ratios 

Before discussing the implications of the observed low 
[a/Fe] ratios in terms of a nucleosynthetic origin, we must 
first rule out other (non- intrinsic) systematics. In Paper 
I, we have already discussed the issue of dust and pho- 
toionization in the double DLA toward Q2314-409. Here, 
we embark upon a similar assessment for the CTQ247 
triplet. In particular, we investigate whether systematics 
may cause (a) the a/Fe ratio to be lower in MDLAs than 
in single DLAs; and (b) the a under-abundance in MDLAs 
to be more evident in S rather than Si. 

First, we consider dust depletion. Although disentan- 
gling dust depletion from pure nucleosynthesic effects in 
DLAs is still matter of debate, some broad trends can 
be established: (1) In no instance is S known to be de- 
pleted in the Galactic ISM (e.g. Savage & Sembach 1996). 
Conversely, Si is mildly depleted in most environments and 
Fe is very easily incorporated into dust. Therefore, dust 
depletion will cause an increase in both [S/Fe] and [Si/Fe], 
the effect being most pronounced in the former ratio, so 
[S/Si] > 0. (2) Zn is undepleted; thus, assuming the same 
nucleosynthetic origin as Fe, one expects [Zn/Fe] > in 
dusty systems. 

Both in CTQ247A and B the [Zn/Fe] shows a signifi- 
cant departure from solar, indicating that dust is present. 
This is also true for Q2314-409, the other MDLA for 
which a Zn abundance is determined; in all cases [Zn/Fe] > 
0.3, typical of other DLAs. Therefore, any downward cor- 
rection due to dust would further emphasize the low a/Fe 
in these MDLAs. It is noteworthy that Q2359-02A, which 
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has the largest [Si/Fe] of all, also shows an unusually large 
dust content, [Zn/Fc]=+0.88, explaining its anomalously 
large [Si/Fe]. Not only is Q2359-02A relatively highly de- 
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Fig. 4. UVES spectrum showing nonsaturated transitions 
in CTQ247C (z abs = 2.6215) 

pleted compared to other MDLAs, but also in comparison 
with other DLAs which mostly have [Zn/Fe] < 0.6. 

Can low dust content in MDLAs compared with sin- 
gle DLAs be responsible for lower [S,Si/Fe] in the for- 
mer? Fig. HJlshows [S/Si] vs. [S/Fe] for the (few) systems 
where both elements are available, including 3 MDLAs. 
It can be seen that apart from one DLA with an ex- 
tremely large [S/Fe], there is no strong trend for sys- 
tems with lower [S/Si] (supposedly less dusty) to show 
lower [S/Fe], suggesting that low S/Fe ratios are not 
strongly effected by dust. Indeed, the large S abundance 
measured in Q0307— 49 might be due to Lya blending 
( Bonifaci o et al. 2001(1 and must be considered with cau- 
tion. Therefore, low S/Fe is not obviously correlated with 
low dust content, at least not to the accuracy we can 
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Fig. 5. Relative ratios of a/Fe for MDLAs (solid circles), 
the DLAs toward Q0201+1120 and Q0000-262 (solid tri- 
angles) and DLAs taken from the literature (stars). All 
ratios have been corrected to the following solar val- 
ues: log (S/H) Q = -4.80, log (Si/H) Q = -4.44, log 
(Fe/H) = -4.50 (Grevesse & Sauval 1998). The data 
points shown in this plot are listed in Table |3I 
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Fig. 6. [S/Si] vs. [S/Fe] in MDLAs (circles), single DLAs 
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presently measure this ratio. In conclusion, it does not 
seem from the present sample that a general trend of low 
a/Fe values be driven by any effect of dust, nor do MDLAs 
have atypical dust content. 
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Table 3. Abundance ratios: ° ]Prochaska et al.~( 2001) 
fc Prochaska k Wolfe (1999); c Lu, Sargent & Barlow (1998); 
d Lu et al. (1996); e Dessauges-Zavadsky et al. (2001); 
^Lopez et al. (1999); 9 Lopez et al. (2002); ^Bonifacio et al. 
(2001); Tevshakov et al. (2002); ^Ge, Bechtold & Kulkarni 
(2001); fc Pettini et al. (2002); 'Prochaska, Gawiser & Wolfe 
(2001); "'Centurion et al. (2000); "Prochaska & Wolfe 
(1996); "Dessauges-Zavadsky, D'Odorico & Prochaska (in 
prep.); p This paper; 9 Paper I 



We now turn to ionization. Given the large A^(Hi) 
of CTQ247A and B, it is highly unlikely that ioniza- 
tion corrections are required l |Viegas 1995| ). Nonetheless, 
Prochaska et al (2002c) found significant corrections in a 
DLA with log A^(H i) = 20.8, so there are occasional cases 
where photoionization plays a part even in relatively large 



column density systems. In order to constrain the ion- 
ization parameter for CTQ247A and B we use CLOUDY 
<|Ferland 1993jl models with a log N(H i) = 21.1, a Haardt- 
Madau ionizing spectrum IjHaardt fc Madau 1996|l at z — 
2.3 and a metallicity of —1.5 solar, and use the ob- 
served ratios to constrain the ionization parameter. For 
CTQ247A, R = logAT(Alm)/iV(Alii) provides the most 
stringent limit on log U, whereas the ratio of Fein/Fen 
provides most information on the ionization of CTQ247B. 
In both cases, the limiting ionization parameter is found to 
be logJ7 < —3.2, see upper panel of Fig. [7J Consequently, 
this limits the correction to [S/Fe] to be less than 0.07 dex. 
Note that a correction this extent is in agreement with the 
subsolar [S/Si] observed for the 3 MDLAs in Fig. and 
brings this ratio into line with solar, within the measure- 
ment errors. 

CTQ247C has a significantly lower Hi column density. 
Our spectrum does cover the Fein A1122 transition but 
the limit is not very meaningful. However, we can com- 
pare Al II and Al in whose ratio has been shown to exhibit 
a steady trend with AT (Hi), a higher fraction of Aim be- 
ing present for lower AT (Hi) iVlad ikTet al. 200T)l . We find 
R = —0.55 which is typical of moderate column density 
DLAs, log AT (Hi) > 20.5. Repeating the above CLOUDY 
models with the appropriate Hi confirms that the cor- 
rections to [S/Fe] will not exceed 0.1 dex in this MDLAs 
(bottom panel in Fig. |7J). Therefore, we would need to 
make a small upward correction to [S/Fe] of up to about 
0.1 dex but this does not alter our conclusions. 

We note that the MDLAs toward Q2359-02 have 
both high R w —0.3, indicating relatively high ionization. 
According to our CLOUDY models, however, the correc- 
tions for [S/Fe] are still negligible for the corresponding 
ionization parameter, while [Si/Fc] requires downward cor- 
rections of 0.1-0.2 dex. 

From a more general point of view, Fig. dtells us that 
for a typical log N (Hi) = 20.5 DLA the corrections are 
downward for Si/Fe and upward for S/Fe, whereas at the 
low end of Hi column densities the [Si/Fe] ratio is signif- 
icantly more sensitive to ionization than [S/Fe] and both 
corrections are downward (cf. Procha ska et al. 2 002c). 
Ionization will therefore require downward corrections in 
[Si/S] but these will not necessarily contribute to the wider 
spread of MDLA [Si/Fe] values we observe in Fig. O 

We can quantify the effect of ionization on the compar- 
ative dispersion between [Si/Fe] and [S/Fe] by performing 
Monte-Carlo simulations with a grid of CLOUDY mod- 
els. To this end, we constructed 2 samples of A^ = 500 
CLOUDY models each, with -4 < logU < -3 ('low- 
ionization' sample) and —4 < \ogU < —2 ('high+low- 
ionization' sample), and Hi column density and metal- 
licity values uniformly distributed in the ranges 19.4 < 
logAT(Hi) < 22.0, -2.4 < log(Z/Z ) < -0.8 5 . To simu- 
late our observations we selected randomly 1 000 sets of 20 

5 Of course, uniformly distributed values of H I and Z do 
not represent the observed distribution of DLAs, this range is 
purely for the purpose of investigating the effects of ionization. 
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Fig. 7. Ratio of doubly to singly ionized species as a 
function of ionization parameter U from CLOUDY models 
with a Haardt & Madau l|Haardt fc Madau 1 996 ) ionizing 
spectrum. Top panel: Model with logTV(Hi) = 21.1 and 
gas metallicity \ogZ = log^Q — 1.5 compared with data 
points for CTQ247A (open square) and CTQ247B (filled 
squares). Middle panel: The measured Alin/Aln ratio 
for CTQ247C, and a model using logiV(Hi) = 20.5 and 
metallicity logZ = logZQ — 2.5. Bottom panel: A model 
using logiV(Hi) = 20.0, the column density cutoff in our 
definition of MDLAs, and metallicity log Z — log Zq — 1.5. 
In all panels the dashed lines indicate solar-corrected ra- 
tios while the solid lines are for column density ratios. 



DLAs each from the two samples. Their [X/Fe] distribu- 
tions are shown in the upper panel of Fig. [HI The difference 
between the scatters in [Si/Fe] and [S/Fe] was calculated 
as AS = a [Si/Fe] - <r[S/Fe], where cr[X/Fe] is the stan- 
dard deviation of each dataset. The bottom panel of the 
figure shows the AS distribution for both samples. For the 
low+high-ionization sample the probability P(AS > 0) 
of having a wider scatter in [Si/Fe] is roughly 0.7, while 
P(AS > 0) = 0.85 if the DLAs are drawn exclusively from 
the low-ionization sample. This is an indication that ion- 
ization does affect the distributions differentially, although 
the differences are generally below measurement errors. 

4.3. Other possible systematic effects 

Finally, other possible systematic effects we have to be- 
ware of are: (1) Column densities: our column densities 
come from x 2 fits whereas the majority of literature DLAs 
have AODM measurements. However, it is well established 
that both approaches give essentially same results pro- 
vided the lines are nonsaturated as is our case. (2) Solar 
abundances: solar abundances have undergone some sig- 
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Fig. 8. Simulated [S/Fe] and [Si/Fe] distributions (top 
panel) and difference between their scatters (bottom 
panel). The solid (dotted) histogram shows DLA simula- 
tions drawn from a sample of DLAs with — 4 < log U < — 3 
(-4<logC7<-2 ) 



nificant revisions in recent years; all the points in figures [5] 
and[!|]have been put on a common solar scale. (3) Atomic 
data: some /-values have changed in recent years, most 
notably Fen A1611. However, this should not be an issue 
as long as very few systems rely on only one transition 
6 . (4) Blending: Sn transitions are normally in the forest 
and therefore particularly susceptible to contamination. 
This is especially a worry if the AODM is used since it 
would cause larger N(S) if contamination is present. Since 
40% of literature S measurements comes from the UCSD 
database which uses the AODM, in theory we could have 
compared our data with overestimations of S/Fe. However, 
checking in Prochaska et al. (2001) on a case-by-case ba- 
sis shows that this is unlikely to be a problem because the 
possibility of contamination was carefully assessed in each 
case. Even excluding [S/Fe] points determined by AODM 



e.g., Q2206— 19A is the only DLA in the sample that relies 
solely on Fell A1611 



10 



Sebastian Lopez and Sara L. Ellison: Abundances in Multiple DLAs 



from Fig. [S] still shows the low S trend in MDLAs. (5) 
Instrument: almost all literature abundances come from 
HIRES data whereas we have used UVES at a slightly 
higher resolution. However, some high a/Fe UVES values 
in Table and the low a/Fe HIRES value in the MDLAs 
toward Q2359— 02 show that a possible bias due to differ- 
ent data type is unfounded. 

In summary, we have found that MDLAs exhibit signif- 
icantly lower [S/Fe], and to a less striking extent, [Si/Fe] 
ratios than single DLAs taken from the literature. We have 
argued that these low values are not driven by ionization 
differences or dust depletion, although an observational 
bias may be partially responsible for the clearer abun- 
dance distinction in S compared to Si. Having excluded 
the main systematic effects, these results therefore imply 
a nucleosynthetic origin for differences in abundance ratios 
between MDLAs and single DLAs. 



5. Discussion 

5.1. The Nature of MDLAs 

If the typically low values of [a/Fe] of MDLAs are 
truly nucleosynthetic in origin, this may indicate that 
these systems either represent a particular galaxy type 
or that they share a common evolutionary or environmen- 
tal link. The first issue to address is whether the 'mul- 
tiplicity' is indeed due to more than one distinct galaxy. 



Schaye (2001) has suggested that some DLAs may be as- 



sociated with winds, in which case MDLAs may be caused 
by outflowing material from the parent DLA. Indeed, 



Nulsen, Barcons & Fabian (1998) have shown that large 



scale galactic winds can produce column densities equal 
to those of DLAs. This could be tested by studying abun- 
dances in systems that are effectively transverse multiple 
systems, that is DLAs observed in more than one line of 
sight separated by many tens of kpc. The closest exam- 
ple we have of this in the current dataset is the case of 
Q0201+1120 IjEllison et al. 2001ajl which is located in a 
concentration with at least four other galaxies and indeed 
shows low [S/Fe], possible evidence that the abundance 
trend is symptomatic of a group environment rather than 
simply line of sight velocities within a single object. 

However, this evidence remains vague in our context 
of MDLAs as these imaging studies deal only with DLAs 
without line-of-sight companions. Deep imaging and sub- 
sequent spectroscopy of the fields toward CTQ247 and 
Q2314-409 as well as control fields around QSOs har- 
boring single DLAs is therefore required in order to 
pursue a comparison of possibly different environments. 
At higher redshifts, finding evidence for cluster envi- 
ronments of single DLAs (and even the DLA absorbers 
themselves) has proven difficult IjOawiser et ah 2 001 
IProchaska et al. 2002af) . so much work is still to be done 
in this direction. Similarly, although several DLA imaging 
studies have been carried out (e.g. JLe Brun et al. 19 97). 
we lack studies of the clustering properties in the DLA 



fields at z < 1.0 and there are currently no known MDLAs 
at low redshift. 



5.2. Truncated star formation in MDLAs? 

As suggested in Paper I, the low a/Fe abundances ob- 
served in MDLAs might be due to suppressed SF in the 
galaxies harboring the absorbing gas. We now raise the 
question as to whether environment may cause these ob- 
served abundance anomalies and speculate upon the pos- 
sibility that MDLAs arise in 'protogalactic groups'. 

There is a clear trend of galaxies in rich environments 
out to z 0.5 to exhibit suppressed SF, not just in dense 
cluster cores, but also in clusters out to 2-3 virial radii 
( Bal ogh et al. 19971 Lewis et al. 2002) and perhaps out to 
several Mpc (Pimbblet et al 2001). The line diagnostics 
of these galaxies indicate that radial trends are consis- 
tent with an age sequence in which the last episode of 
star formation happened most recently in the outermost 
galaxies (Ba logh et al. 19991 ). However, the physical mech- 
anisms governing this effect are not clear; violent processes 
such as ram pressure stripping (e.g., ICouch et al. 200ll 
and references therein) and passive exhaustion of gas sup- 
ply (e.g., Balogh et al. 1999| ) are the two main possibil- 
ities. The recent finding by Lewis et al. (2002) that lo- 
cal density is the dominant parameter in suppressing SF 
rates indicates that extreme cluster-scale processes such 
as ram-stripping by the ICM play a minor role. These 
results also point to the possibility that reduction in 
SF need not take place within the gravitational bounds 
of existing virialized clusters and may occur already in 
looser groups before they are accreted onto larger struc- 
tures. This ties in with cluster histories at lower redshifts, 



because although Dressier et al. (1999) find a significant 
post-starburst population in clusters, there is evidence 
that this fraction is not greater than among the field pop- 
ulation dBalogh et al. 1999| ). 

If the same physical processes are on-going at high 
redshifts, then suppressed SF may be evident in galaxies 
associated with early proto-clusters, such as those identi- 
fied by Steidel et al (1998), even though canonical, viri- 
alized clusters have yet to form. Although the abundance 
ratios for MDLAs that we have discussed here are distinct 
from 'field' DLAs, they have typical metallicities for their 
redshift (c.f. Pettini et al 1999) and moderate dust con- 
tent. These properties are indications that some chemical 
maturity has indeed taken place, and MDLAs are simi- 
lar to other DLAs in terms of their general enrichment. 
However, the low [a/Fe] is an indication that what sets 
MDLAs apart is that they have been evolving quiescently, 
without any further major SF episodes, over the last ~ 
500 Myrs. 

5.3. Other evidence for quiescent star formation 

If the low a/Fe ratios observed in MDLAs are to be 
explained by quiescent SF, then we should expect to 
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Fig. 9. The N/a ratio in MDLAs (circles), 'field' DLAs 
(stars) and Hn regions (dots; see Pettini et al. 2002 for 
references). The triangles are for the 2 DLAs with com- 
panion LBGs. 



see its signature in other elemental ratios that also 
act as 'cosmic clocks'. The most promising possibil- 
ity is the N/a ratio which traces the primary and 
secondary contributions from stars of different masses 
(Henry, Edm unds and Koeppen 2000| ). Fig. shows the 
N/a ratio in MDLAs, 'field' DLAs and H II regions (Pettini 
et al. 2002). We only used measurements of O or S in the 
plot and corrected S values by the solar S/O ratio (1.54 
dex). Prochaska et al. (2002b) have suggested a bimodal 
distribution of N/a and although this needs to be con- 
firmed by more data, it is clear that a number of DLAs 
exhibit low ratios. For MDLAs, the ratios are all high, 
with no points (N/O) < —1.9, supporting the quiescent 
SF scenario. It perhaps indicates more chemical maturity 
than average, but also that the role of massive stars in 
polluting the ISM of MDLAs may be less important. 

The other established trend with metallicity in 
Galactic stars is the odd-even effect. Although this 
is not technically a chronometer in the same sense 
as N and the alpha elements, there is a correlation 
with metallicity which may be due to metallicity de- 
pendent yields IjGoswami fe Prantzos 20 00). At our dis- 
posal we have [Mn/Fe] = -0.12 ± 0.08 (CTQ247A) 
and [Mn/Fe] = -0.31 ± 0.03 (CTQ247B), the for- 
mer of which is among the highest values in DLAs 
of this metallicity ( |Ledoux, Bergeron fe Petitjean 2 002 



Dessauges-Zavadsky et al 2002). If these high ratios are 
not due to dust depletion (Fe is more depleted in the local 
ISM than Mn) or ionization, they show a mild odd-even 
effect. Analogously, [Al/Si] = 0.07±0.12 in CTQ247C, the 
highest value ever measured, and [Al/Si] = —0.28 ± 0.18 
in B2314— 409B (Paper I), which is among the highest 
values after the z = 2. 154 MDLA toward Q2359-02 
l|Prochaska fe Wolfe 2002jl . Altogether, the odd -even ef- 
fect seems to be mild in those MDLAs where we can 
measure it. Since models of massive star yields give good 
agreement with observations of Mn/Fe in Galactic stars 
IjGoswami fe Prantzos 20001 and references therein) and 
also DLAs ( |Ledoux, Bergeron fe Petitjean 2002| ), the lack 
of odd-even effect at low metallicity in MDLAs indicates 
again that high-mass stars have not been the dominant 
metal pollutant in MDLAs. This may be interpreted as 
further evidence that there have not been several epochs 
of SF building on the enriched interstellar media from pre- 
vious episodes; further support of suppressed SF. 

6. Summary 

We have studied a sample of 7 DLAs with line-of-sight 
companions, 'multiple' DLAs in our nomenclature, and 2 
DLAs arising in transverse groups. We have shown that 
the relative abundances in these two sub-samples are un- 
usual as compared with single DLAs. In particular, the 
[S/Fe] and [Si/Fe] ratios are statistically lower than liter- 
ature DLAs, although the effect is more striking in [S/Fe]. 
We suggest that this difference is due to (1) an obser- 
vational bias against measuring low S column densities 
because weak Sn are usually lost in the forest; and (2), 
to a lesser degree, the larger scatter induced in [Si/Fe] 
by photoionization. Besides low a/Fe ratios, MDLAs also 
exhibit a mild odd-even effect and relatively high N/a 
ratios. We interpret these results as truncated star forma- 
tion in MDLAs. If the multiplicity of these DLAs is due 
to grouping in physical space, a thesis we cautiously sup- 
port, environment may be the cause of the quiescent SF 
scenario, just as is observed in more nearby clusters and 
groups of galaxies. 
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